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We discuss the features of instabilities in asymmetric nuclear matter, in particular 
the relation between the nature of fluctuations, the types of instabilities and the 
properties of the interaction. We show a chemical instability appears as an in- 
stability against isoscalar-like fluctuations. Then starting from phenomenological 
hadronic field theory (QHD), including exchange terms, we discuss the symmetry 
energy and the relation to the dynamical response inside the spinodal region. 

1 Introduction 

Nuclear matter (NM) is a binary system (neutrons and protons), this feature 
has stimulated a quite exciting research field: nuclei far from /3-stability, dy- 
namical effects of a large charge asymmetry, properties of neutron stars. In 
particular with radioactive beams (not only), NM with a high isospin asymme- 
try can be created transiently, offering the possibility to study also chemical 
instability associated with isospin asymmetry of NM. Indeed, recent experi- 
ments on nuclear reactions involving different ratios of the neutron to proton 
numbers have shown that products from nuclear multifragjmcntation depends 
strongly on the isospin asymmetry of the colliding nucleian. 

In the first part of the paper we will focus on the behaviour of Asym- 
metric NM (ANM) at sub-normal densities where we can foresee scenarios for 
a dynamical formation of fragments with particular isotopic contents due to 
instabilities. To study these problems we show, by linking thermodynamics 
and Fermi liquid theory, the relation between the nature of fluctuations and 
the types of instabilities for a general binary system. In the second part we 
show a relativistic field model of hadrons (QHD) based on strongly interacting 
nucleon and meson fields. Inside the QHD model we have recently proposed a 
method to introduce Fock term contributions 0. Here we want to stress some 
implications of that improvement on the isospin physics discussed in the first 
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part. Comparisons with the QHD model in Hartree approximation and with 
Skyrme interactions will be addressed. 



2 Chemical and mechanical instabilities 

2.1 Analytical analysis 

In the context of multifragmentation the instability of NM (region where the 
system becomes unstable against long wave length but small amplitude fluc- 
tuatiaK^iilaiys a crucial role. These aspects were discussed repeatedly in the 
past,u'BBLru but the binary character will induce new features for both scenar- 
ios that are absent in one-component systems. In the framework of Landau 
theory for two component Fermi liquids the spinodal border was determined 
by studying the stability of collective modes described by two coupled Landau- 
Vlasov equations for protons and neutrons. In terms of the appropriate Landau 
parameters the stability condition can be expressed asS, 

(1 + F "")(l + F PP ) - F^ p Fq 11 > . (1) 

It is possible to showj-that this condition is equivalent to the following thermo- 
dynamical condition Ej 



dP\ fdfip 
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> . (2) 



discussed intLSII3, where y is the proton fraction. In Fig. 1 we show the spinodal 
lines obtained from eq. (Q) (continuous line with dots) which for asymmetric 
nuclear matter is seen to contain the lines corresponding to " mechanical insta- 
bility", (^Qf^j < (crosses). Therefore both eqs. (|j]]|) describe the whole 

region of instability. We want to stress, however, that by just looking at the 
above stability condition we cannot determine the nature of the fluctuations 
against which a binary system becomes chemically unstable. Indeed, the ther- 
modynamical condition in eq. (^) cannot distinguish between two very different 
situations which can be encountered in nature: an attractive interaction be- 
tween the two components of the mixture {F^^Fq 71 < 0), as is the case of 
nuclear matter, or a repulsive interaction between the two species. 

We define as isoscalar-like density fluctuations the case when proton and 
neutron Fermi spheres (or equivalently the proton and neutron densities) fluc- 
tuate in phase and as isovector-like density fluctuations when the two Fermi 
sphere fluctuate out of phase. Then it is possible to provelij, based on a ther- 
modynamical approach of asymmetric Fermi liquid mixtures, that chemical in- 
stabilities are triggered by isoscalar fluctuations in the first, i.e. attractive, 
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Figure 1: Spinodal lines corresponding to chemical (circles) and mechanical (crosses) insta- 
bility for three value of proton fraction y. 



situation and by isovector fluctuations in the second one. For the asymmet- 
ric nuclear matter case because of the attractive interaction between protons 
and neutrons the phase transition is thus due to isoscalar fluctuations that in- 
duce chemical instabilities while the system is never unstable against isovector 
fluctuations. Of course the same attractive interaction is also at the origin 
of phase transitions in symmetric nuclear matter. However, in the asymmet- 
ric case isoscalar fluctuations lead to a more symmetric high density phase 
everywhere under the instability line defined by eq. (m u 





Figure 2: Time evolution of the density p{x, y) in the plane 2=0 for initial density p(°) = 
0.09/m — 3 , at T=5 MeV and asymmetry / = 0.5. Upper panels show contour plots of p(x,y) 
and lower panels the corresponding two dimensional surface. 
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2.2 Numerical results: heated nuclear matter in a box 



The previous discussion is based on a link between thermodynamics and Fermi 
liquid approaches. Then numerical approach were performed in order to follow 
all stages of fluctuation development in the fragment formation process. In 
the numerical approach we consider nuclear matter in a box of size L=24fm 
imposing periodic boundary conditions. 

We follow a phase space test particle method to solve the Landau- Vlasov 
dynamics (using Gaussian wave packets) so the dynamics of nucleon-nucleon 
collision is also included. An initial temperature is introduced by distributing 
the test particle momenta according to Fermi distribution. We have followed 
the space-time evolution of test particles for a value of the initial asymmetry 
1=0.5, at initial density = 0.09/to~ 3 and T=5 MeV in such way we start 
from the region of chemical instability (see Fig.l). The initial density per- 
turbation was created automatically due to the random choice of test particle 
positions. We report in Fig. 2 density distribution in the plane z—0 at three 
time steps t=0, 100,200 fm/c, corresponding respectively to initial conditions, 
intermediate and final stages of the spinodal decomposition (SD). The contour 
plots delimit the region with density higher than the initial value of density. 
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Figure 3: Time evolution of neutron (thicks lines) and proton (thin lines) abundance (a) and 
asymmetry (b) as function of density. 

We report the time evolution of neutron (thick histogram) and proton 
(thin histogram) abundance (Fig. 3a) and asymmetry in various density bins 
(Fig. 3b). The dashed line respectively shows the initial uniform density value 
p(°) = 0.09/to~ 3 (3a) and the initial asymmetry 1=0.5 (3b). The drive to 
higher density regions is clearly different for neutrons and protons: at the end 
of the dynamical clustering mechanism we have different asymmetries in the 
liquid and gas phases (see the panel at 200 fm/c). This result is the same 
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of ref.10 where the dynamics of mechanical instabilities were studied, demon- 
strating that also in a complete dynamical calculation the kind of fluctuation 
associated with chemical and mechanical instabilities are the same. 

3 Quantum hadrodynamics model (QHD) 

3. 1 Treatment of Fock terms 

Phenomenological hadronic field theories (Quantum Hadrodynamics, QHD) 
are widely used in dense nuclear matter studies E3. In most of the previous 
works on the subject, the Relativistic Mean Field (RMF) approximation of 
QHD has been followed. In the RMF the meson fields are treated as classi- 
cal fields and a Hartree reduction of one body density matrices is used. This 
implies that each meson field is introduced, with appropriated readjusted cou- 
plings, just to describe the dynamics of a corresponding degree of freedom, 
without mixing due to many-body effects: neutral a and lo mesons are in 
charge of saturation properties, isospin effects are carried by isovector 5 and 
p mesons. In a sense the model represents a straightforward extension of the 
One-Boson-Exchange (OBE) description of nucleon-nucleon scattering. 

Our aim is to introduce explicit many-body effects just evaluating exchange 
term contributions. Fock terms play an essential role in symmetry breaking 
and consequent mixing of different degrees of freedom. In particular, in the 
context of the QHD model, essential properties of nuclear matter come mostly 
from the two neutral strong meson fields. Hence it is important to evaluate 
the Fock contribution associated with these fields. ._. 

Then we will start from a QHD — II model tL3 where the nucleons are 
coupled to neutral scalar a and vector u mesons and to the isovector p meson. 
Sclf-intcractionj-tcpps of the u-ficld were originally introduced for renormal- 
ization reasons LLjEJj and can also be considered as a way to parametrize the 
density dependence of NN force. Actually they are also describing medium ef- 
fects essential to reproduce important properties (compressibility and nucleon 
effective mass) of nuclear matter around saturation density. 

The Lagrangian density for this model is given by: 

C = V47„(^ - 0yV - g p B» ■ t) - (M - 5s 0)]^ + ^(d^d^ - m s V) 

-^0 3 - \^ - \w» v W^ + \m v 2 V v V v - ' L ^ + \™ P 2 ®» ■ B" (3) 

where W» v (x) = d»V v {x)-d v V»(x) and L^(ar) = d^B v (x)-d v B> i {x) . Here 
-0(x) generally denotes the fermionic field, </>(x) and V(x) represent neutral 
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scalar and a vector boson fields, respectively. B v {x) is the charged vector field 
and t denotes the isospin matrices. 

We have treated Hartree-Fock (HF) terms in EOS and transport equation 
at the same level, for this reason we performed the many-body calculations in 
the quantum phase space introducing the Wigner transform of the one-body 
density matrix of the fermion field. This method has the advantages of a direct 
derivation of dynamical transport equations 113. In ref. u the way we introduce 
Fock terms is discussed in detail, here we want to discuss some results for EOS 
of ANM. 




Figure 4: Symmetry energy per nucleon vs. baryon density. Long dashed line: NLH with p 
meson. Solid line: NLHF results. Dashed line: NLH with p and 8 mesons. 



3.2 Asymmetric Nuclear Matter 

We obtain, thanks to Fock terms, scalar and vector isovector contributions to 
symmetry energy, generally associated respectively with 5 and p mesons, even 
without isovector mesonstSI. We show the comparison between our Non Linear 
Hartree-Fock (NLHF) and those of the Non linear Hartree (NLH), including 
the isovector p and S mesons Ej, with parameters fitted in order to give the 
same saturation properties and a symmetry energy at po of 31 MeV. The region 
in (T,p) plaii^at different asymmetry of chemical and mechanical region are 
quite similarEy to the one of Skyrme interaction (see Fig.l). 

In Fig. 4 the symmetry energy obtained in NLHF (solid line) in comparison 
with the result of NLH including both p and S meson (dashed line) . Actually 
the kinetic contribution is subtracted (anyway it is the same in all the mod- 
els), so we will refer to potential symmetry energy EP°!^. For reference the most 
common result, including only p as isovector meson, is plotted (long dashed 
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line). In all these relativistic models a quite repulsive density dependence of 
the simmetry term is obtained, but we notice that the density dependence of 
Efym m t ne complete NLH+p+<5 model is quite different respect to our result. 
This is due to the fact that in NLHF the coupling in the isovector channels 
become density dependent. This implies a "softer" behaviour of the potential 
symmetry term below the saturation density in the NLHF case and a "stiff" 
behaviour more similar to NLH +p above. 




0.0 0.2 0.4 0.6 0.8 1.0 
I=(N-Z)/A 



Figure 5: Ratio of isovector and isoscalar amplitudes as function of asymmetry I at half 
saturation density. The lines have the same meaning of Fig. 4 



A transport equation can be consistently derived to be used for the study 
of dynamical evolution of nuclear matter far from normal conditions. We 
expect that the dynamical evolution can point out the important difference 
between NLH and NLHF by comparison with experimental data for isospin of 
IMF (in multifragmentation or neck events) and for isospin flows. We have 
just found as the latter are affected by density behaviour of E v s '^ rn and the 
nucleon effective mass splitting c3. In this view ve have performed the study 
of collective response in ANM by means of a relativistic kinetic equation in 
linear approximation respect to the fluctuation E3. In Fig. 5 we present the 
unstable isoscalar-like solution of the dispersion relation for low density NM 
(the density is p = 0.4po) 7 in particular the ratio Sps/Sp^ as function of the 
initial asymmetry. We notice that the NLHF tendency to restore the isotopical 
symmetry in the spinodal decomposition of neutron rich NM is different respect 
to NLH result (both dashed and long dashed result). In particular we stress 
that the discrepancy is equal to the one obtained E3 with Skyrme interaction, 
in a non relativistic approach, comparing a Asy-soft with an Asy-stiff term 
for symmetry energy. However at variance with Skyrme results the behaviour 
of E sym in NLHF appear to be Asy-soft at subnuclear density but it's Asy- 
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stiff above saturation density. A mixed behaviour is predicted with expected 
interesting effects on experimental observable depending on the probed barion 
density region of the interacting nuclear matter. 
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